Spin-Isospin Structure and Pion Condensation in Nucleon Matter 
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Abstract 



We report variational calculations of symmetric nuclear matter and pure 
neutron matter, using the new Argonne wis two-nucleon and Urbana IX three- 



• nucleon interactions. At the equilibrium density of 0.16 fm the two-nucleon 

^ . 

. densities in symmetric nuclear matter exhibit a short-range spin-isospin struc- 

!>■ ■ 

ture similar to that found in light nuclei. We also find that both symmetric 
nuclear matter and pure neutron matter undergo transitions to phases with 

, 

^ . pion condensation at densities of 0.32 fm ^ and 0.2 fm ^, respectively. Nei- 

ther transtion occurs with the Urbana V14 two-nucleon interaction, while only 
the transition in neutron matter occurs with the Argonne vn two-nucleon in- 
teraction. The three-nucleon interaction is required for the transition to occur 
in symmetric nuclear matter, whereas the transition in pure neutron matter 
occurs even in its absence. The behavior of the isovector spin-longitudinal 
response and the pion excess in the vicinity of the transition, and the model 

dependence of the transition are discussed. 

PACS numbers: 21.65.+f, 21.30.Fe, 26.60.-hc, 64.70.-p 
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I. INTRODUCTION 



A central problem in nuclear many-body theory is the prediction from realistic models 
of nuclear forces, of properties of uniform nuclear matter in stable nuclei and in neutron 
stars. This problem has been the focus of much attention since the pioneering work of 
Brueckner, Levinson and Mahmoud in 1954. Early work in this field, based primarily 
on Brueckner-Bethe-Goldstone many-body theory, was reviewed by Bethe in 1971. A 
variational method based on hypernetted chain summation techniques (VCS, previously 
denoted FHNC-SOC) was developed in the 1970s [0, particularly to include the effects of 
many-body correlations, presumably important in dense neutron star matter. Calculations 
performed since then have confirmed that many-body clusters make significant contributions 
to the binding energies of equilibrium nuclear matter and light nuclei 

Friedman and Pandharipande (hereafter denoted FP) carried out detailed calcula- 
tions of the equation of state of symmetric nuclear matter (SNM) containing equal numbers 
of protons and neutrons and of pure neutron matter (PNM) with VCS, using the Urbana 
fi4 model of the NN interaction [|10|. In these calculations, effects of many-body forces 
were introduced partly via a small density dependence in the NN interaction and partly 
by adding a density-dependent contribution, attributed to the attractive two-pion exchange 



three-nucleon interaction W^. The parameters of these contributions were adjusted to repro- 
duce the empirical properties of SNM, and methods were developed to interpolate between 



SNM and PNM ||12|. The main weakness in this approach is that the effect of many-body 
interactions on the structure of matter is ignored. 

In the 1980s, a variational theory using Monte Carlo methods (VMC) was developed 



for light nuclei [|l^. Correlation operators with the same structure were employed in both 
VMC and VCS methods, with the aim of obtaining a unified description of light nuclei and 
nuclear matter. The Urbana VII (UVII) model of three nucleon interactions, containing the 



Fujita-Miyazawa two-pion exchange three-nucleon interaction |TJ] and a phenomenological, 
shorter range interaction, was added to the Hamiltonian. The parameters of the UVII were 
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determined by reproducing the binding energies of ^H, ^He and the equihbrium density of 
SNM using approximate VMC and VCS calculations |T5[. 

Wiringa, Fiks and Fabrocini W^ (hereafter denoted WFF) calculated the equation of 



state of SNM and PNM with an improved VCS method using the Urbana and Argonne ||T7 
f 14 models of the NN interaction and the Urbana VII three-nucleon interaction. Their calcu- 
lations of PNM indicated a transition to a new phase, possibly including pion condensation, 
at p ~ 0.24 fm~^ with the Argonne fi4, though not with Urbana ^14. No evidence of a 
similar transition in SNM was found by WFF. 

Since the pioneering work of Migdal 0] and of Sawyer and Scallapino |]l9l , many investi- 
gators have used effective interactions to study the possibility of pion condensation in SNM 



and PNM. These efforts were recently reviewed by Kunihiro et al. [^. In these studies, the 
A-resonance is explicitly considered as a non-nucleonic degree of freedom. In contrast, WFF 
considered only nucleonic degrees of freedom and absorbed the effect of the A-resonance into 
the two- and three-nucleon interactions. The approaches used by Migdal and by WFF to 
study pion condensation are thus very different. Although the effect of this type of transi- 
tion on the equation of state is relatively small, it can have important consequences for the 
cooling and evolution of neutron stars . 

In the late 1980s Carlson [H developed the Green's function Monte Carlo (GFMC) 
method, with which exact calculations of light nuclei are possible, starting from Hamiltonians 
with realistic two- and three-nucleon interactions. The GFMC method can be used to 
determine some of the parameters of the three-nucleon interaction from exact calculations 
of nuclear binding energies. 



Recently, Wiringa, Stoks and Schiavilla made significant improvements to the 
Urbana-Argonne fi4 models of the NN interaction by including isospin symmetry break- 
ing terms. Their resulting Argonne fig is one of the new models fit to the Nijmegen NN 



scattering database p^j25[] . Using this interaction, along with the GFMC method, Pudliner 



et al. obtained parameters for the Urbana IX (UIX) model of the three-nucleon inter- 



action. Exact GFMC calculations by Pudliner et al. [pTtPSj have shown that the Argonne 



Vis and Urbana IX interactions provide a good description of all bound states, up to seven 
nucleons. While all observed energies are not exactly reproduced, in each case the difference 
between theory and experiment is much smaller than the contribution of the three-nucleon 
interaction. 

In the present work, we employ the Argonne Vis and Urbana IX interactions to study 
the properties of SNM and PNM, using the VCS method. We have introduced further im- 
provements to the method, in order to address the relatively strong momentum dependence 
of the Argonne vis interaction. These improvements are described in Appendices A and B. 

With these interactions, we find that both SNM and PNM undergo transitions to a new 
phase at densities of ~ 0.32 fm~^ and ~ 0.2 fm~^, respectively. The transitions are similar 
in nature to that found by WFF in PNM, using the Argonne vu and Urbana model VII. 

In this paper, we discuss the phase transition and the spin-isospin structure of normal 
SNM at equihbrium density. Additional results for the equation of state of dense matter, 
including relativistic effects, will be reported separately. The plan of the present paper is 
as follows. The Hamiltonian and the VCS method are reviewed in Sections II and III. In 
Section IV, we present various pair distribution functions of interest, at both equilibrium 
and transition densities. These results clearly indicate that in normal SNM at p = 0.16 fm~^, 
the tensor correlations in T=0 states have near maximum possible strength at r ~ 1 fm. 
The results also indicate the important role played by the short-range parts of (Tj • (TjTi ■ Tj 
and (Tj • (Tj interactions in determining the transition density, as predicted by Migdal. In 
Migdal's approach, the strength of the short-range part of the CTj • cr^Tj • Tj interaction is 
represented by the Landau- Migdal parameter g'. We also comment on the model dependence 
of this transition in Section IV, by comparing the Urbana vu and the Argonne vu and 
Vis interactions. Results for the sums of isovector spin-longditudinal responses of SNM 
and PNM are presented in Section V. The corresponding results for the strengths of pion 
exchange interactions and pion fields appear in section VI. These results clearly indicate 
that the high density phase includes pion condensation, as anticipated by WFF. Results for 
SNM are presented at normal density and at the transition density. The response of SNM 
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at normal density exhibits little indication of the phase transition, which is evidently of first 
order in the present approach. Conclusions are presented in Section VII. 



II. THE HAMILTONIAN 

The Argonne vis two-nucleon interaction has the form: 

Vl8,ij = + Vem- (2.1) 

p=l,18 

The electromagnetic part Vem consists of Coulomb and magnetic interactions in the nn, np 
and pp pairs, and it is omitted from all nuclear matter studies. The strong interaction part 
of the potential includes fourteen isoscalar terms with operators: 



1, (Ti ■ aj, Sij, (L • S)i„ L\ L'cr, ■ o",-, (L • S)?- ® [1, r, • r,] . (2.2) 



By convention, operators with even p < 14 include the ■ Tj factor, while those with odd 
p < 13 are independent of isospin. The operators with p =15, 16 and 17 are isotensor: 

Ofr^^'^'^ = {"iT^iTzj - Ti ■ Tj) (g) (1, (Ti ■ (Tj, Sij). (2.3) 

The interactions associated with operators p = 16 and 17 are generated by the mass differ- 
ence between neutral and charged pions, while the phenomenological interaction with the 
p —15 operator is necessary to fit the difference between np and pp, ^Sq phase shifts. Finally, 
the isovector interaction with the operator, 

Ojf = (r,, + T,,), (2.4) 

is required to fit the nn scattering length and effective range along with the np and pp data. 
The interaction constructed with the above elements provides an accurate fit to the NN 
scattering data in the Nijmegen database, using the kinetic energy operator: 



( 1 1 \ [ 1 1 , _2 



rUp rrinj \mp rrir, 
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which takes into account the mass difference between protons and neutrons 



V.^ (2.5) 



The Urbana models of Vijk contain two isoscalar terms: 

V^,k = V^l + V,f, . (2.6) 
The first term represents the Fujita-Miyazawa two-pion exchange interaction: 

V^l = J2cyc^27T ({t i ■ Tj,Ti ■ Tk} {Xij, Xjfc} + ^ [t j ■ Tj , Tj ■ T fc] [Xjj , Xjfc]) , (2.7) 

Xij = SijT^{rij) + (Ti ■ ajY^{rij), (2.8) 

with strength A2tt- The functions TT,{rij) and l^,r('^ij) describe the radial shapes of the 
one-pion exchange tensor and Yukawa potentials. These functions are calculated using the 
average value of the pion mass and include the short-range cutoffs used in the Argonne Vig 
NN interaction. The term denoted by V^f^^ is purely phenomenological, and has the form: 

V,% = U,Y.Tl{n^)Tl{nk). (2.9) 

eye 

This term is meant to represent the modification of NA- and AA-contributions in the two- 
body interaction by other particles in the medium, and also accounts for relativistic effects. 
The spin-isospin dependence of these effects is neglected. 

The two parameters A2tt and Uq are chosen to yield the observed energy of and 
the equilibrium density of nuclear matter, po = 0-16 fm^^. The parameters of Model VII, 
used by WFF, were determined from approximate variational calculations using the older 
Argonne Vi^ NN interaction, and have the values: A2t^ = —0.0333 MeV and Uq = 0.0038 
MeV. The parameters of Model IX, = -0.0293 MeV and Uq = 0.0048 MeV, have been 
determined from exact GFMC calculations of and the present approximate variational 
calculations of nuclear matter with the Argonne fig NN interaction. Most of the difference 
between the parameters of Models VII and IX is due to the use of exact, rather than 
approximate calculations of in the latter. The parameters of Model VIII, also determined 
via exact calculations, albeit with the older Argonne fi4 interaction, have the values: 
= —0.028 MeV and Uq = 0.005 MeV, which are not very different from the present 
Model IX parameters. The insensitivity of these parameters to the details of the Urbana- 



Argonne models of Vij had been previously noted p9[. For example, calculations of H and 



nuclear matter saturation density using Model VII yielded similar results with the Urbana 
Vi4 and Argonne fi4 interactions, despite significant differences in the tensor components 
of these interactions. We also note that the values of A2tt obtained from all these fits are 
close to ~ —0.03 MeV, the strength predicted by the Fujita-Miyazawa model. More recent 
theoretical models of V^^l, such as the Tucson- Melbourne predict significantly larger 
strengths, A2Tr ~ —0.063 MeV. These models of also have additional terms which are 
neglected here. 

The isotensor and isovector parts of Vij, and the isovector part of the kinetic energy K, 
are very weak and may presumably be treated as first order perturbations. In first order, 
these terms do not contribute to the energy of SNM, which has total isospin T=0. Therefore, 
the SNM calculations, which neglect Vem by definition, are carried out with the isoscalar part 
of the Hamiltonian: 

HsNM = E -f f + v.^ + E E -'ir^.)of, + E v^.^- (2-10) 

■ ^ yiILp IILnJ i<j p=l,U i<j<k 

The isotensor and isovector terms do contribute to the energy of PNM. The isospin operators 
Tj ■ Tj, Sr^iTzj — Tj ■ Tj and r^j + t^j reduce to constants, respectively equal to 1,2 and —2, 
in PNM. Hence, the full Hamiltonian, excluding Vgm, has a simple form for PNM, given by: 

HPNM = E -^v? + E E v''-'ir.,)o^r + E V^J^> (2-11) 



2"^" j<j g=i J i<j<k 



with 



= t,i + t,2 + 2f - 2t;^^ (2.12) 
= + + 2v^^ (2.13) 
= + + 2v^'' (2.14) 

^2g-l _ ^2g-l ^ ^2g fo^ g > 4 (2.15) 

We note that [xj ■ tj, ti ■ t^] = in PNM, thus the commutator term of V^^l can be omitted 
from YipNM- 
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III. VARIATIONAL CALCULATIONS 



The variational method developed in References [^,11 ,16,3T| provides a means to calculate 
the energy and wavefunction of the ground state of nuclear matter from realistic models 
of nuclear forces. We outline the method here and describe additional developments in 
Appendices A and B. 

The variational wavefunction has the form: 

^. = (^n^^.)'^' (3-1) 

consisting of a symmetrized product of pair correlation operators Fij operating on the Fermi 
gas wavefunction $. In SNM, the F^j include eight terms: 

= E nn^W^,, (3.2) 

p=l,8 

representing central, CTj ■ cr^, Sjj and (L ■ S)jj correlations with and without Tj ■ Tj factors. 
In PNM, the Fij reduce to a sum of four terms with only odd p <7. 

This wavefunction is clearly too simple to accurately describe the ground state of nuclear 
matter. Monte Carlo studies of few-body nuclei [^] and ^^0 0] use additional three-body 
correlation operators in the variational wavefunction. These additions to the wavefunction 
lower the ground state energy of ^^O by ~1 MeV/nucleon. Attempts to include three-body 



correlation operators in variational calculations of SNM are currently in progress |]33|. Exact 
GFMC calculations have now been carried out for nuclei containing up to seven nucleons 
26| , |28[| . The ground state energy of ''Li obtained with variational wavefunctions including 



three-body correlation operators is greater than the exact value by ~0.7 MeV/nucleon. 
From these results, we estimate that the present may underbind SNM by a few MeV. 
In contrast, the three-body correlations have a much smaller effect on the energy of pure 
neutron drops [§. The variational energy of a drop with eight neutrons, calculated with 
the simple is greater than the exact value by ~0.5 MeV/nucleon. Thus, we believe the 
energies calculated in the present work to be relatively more accurate for PNM than for 



8 



SNM. This result is to be expected, since SNM has strong tensor correlations in two-body 
np states, with isospin T=0, and in three-body nnp and pnn clusters with T=i. 

Despite the aforementioned shortcomings, the simple having only pair correlation 
operators describes the gross features of the nuclear wavefunction rather well. For example, 
the spin-isospin dependent two-nucleon distribution functions calculated in this approxima- 
tion are close to the exact distribution functions We note that the magnitude of the 



ground state energy is generally much smaller than those of the positive kinetic and negative 
interaction energies. The magnitude of the error in the ground state energy is typically only 
a few percent of the interaction energy. We can therefore study the main features of the 
spin-isospin structure of PNM and SNM using the simple described above. 



The correlation operators Fij are determined from Euler-Lagrange equations [jTT| that 
minimize the two-body cluster contribution of an interaction {v — X), where: 



Y: a^v^{r,Mv (3-3) 

p=l,14 

A., = E A^(^«.)Of,- (3.4) 

p=l,8 

The variational parameters are meant to simulate the quenching of the spin-isospin 
interaction between particles i and j, due to flipping of the spin and/or isospin of particle i 
or j via interaction with other particles in matter ||3^. We use: 

= 1 for p = 1 and 9 (3-5) 
= a otherwise, (3.6) 

since the operators 1 and (p = 1 and 9) are independent of spin-isospin. The of 
the (L ■ S)]j interactions (p = 13, 14) are also set to unity. The (L ■ S)^^ operator has a 
significant part, which should not be quenched, and quenching the remainder of v^^^^, 
along with v^^^^, does not lower the variational energy. The A^(r) simulate screening effects 
in matter and are determined by the healing distances of the correlation functions f '^: 

Fir > dn = (3.7) 



The radial dependence of A^(r > d^) is determined from the above constraint. The A^(r < (i^) 
are constants chosen to make the gradients of at equal to zero. The energy of nuclear 
matter is particularly sensitive to the range of the tensor and central correlations. For 
simplicity we assume: 

dP = dt for p = 5, 6 (3.8) 
dP = dc for p ^ 5, 6. (3.9) 

The Fij, and consequently "^y, thus depend on three variational parameters: a, dc and dt, 
the values of which are determined by minimizing the energy. Two additional parameters 
used by WFF have only a small influence on the energy. 

It is convenient to divide the interaction and correlation operators as follows: 

% = Vs,ij + Vb^ij + Vq^ij, (3.10) 
-P^ij — Ps,ij ~\~ -^b,ij- (3-11) 

The static parts Vs and Fg involve terms with momentum independent operators Of~^~^. 
The spin-orbit terms, with p =7 and 8 are included in Vf, and Ff,, while terms quadratic 
in L, with p =9 to 14 are included in Vq. The two-body (2B) cluster contributions involve 
expectation values of the operators: 

r(r^,)0'!^v!pf^r'ir,,)0'!;, (3.12) 
^(r^,)0!^Vy'^'(r,J)0<, (3.13) 

in plane-wave states. We divide these contributions into five parts: the static parts, (f )-2B-s 
((T)-2B-s) include interaction (kinetic) energy contributions for q,p and q' {q and q') < 6; 
the spin-orbit parts, {v)-2B-h and (T)-2B-b include terms with p < 8, and one or more of 
the p, q, and q' equal to 7 or 8; and the quadratic part, (f )-2B-q includes all terms with p > 
9. 

The results obtained with various Hamiltonians for SNM at kj=1.33 fm~'^ and 1.6 fm~^ 
are listed in table |. These Hamiltonians include the kinetic energy plus Urbana vi^ vij 
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(U14), Argonne vis Vij (A18) and Argonne vis Vij plus Model IX Vijk (A18+IX). The Fermi 
kinetic energy is listed as (T)-IB, the one-body cluster contribution to the kinetic energy. 
Variational cluster Monte Carlo studies of ®n-drops and ^^O including up to 



five-body cluster contributions to the ground state energy, indicate that the convergence 
of the cluster expansion is not particularly good in either case. For example, the one- 
through five-body clusters contribute 12.9, —54.5, 11.1, —3.8, and 1.1 MeV respectively to 
the ground state energy of the ^n-drop. This result does not contradict the earlier lowest- 
order constrained variational calculations of neutron gas |36|j37|| , which retained only the 



(T)-IB and all 2B contributions. The aim of the earlier studies was to obtain estimates 
of the ground state energy; the dc and df were not determined variationally, but rather 
were fixed by constraints. In the present work, is determined variationally to study the 
structure of SNM and PNM. The larger optimum values of dc and dt lead to the significant 
size of many-body (MB) clusters of three or more particles. 

According to Table |, the kinetic energy and static interactions and correlations make 
large contributions to the IB and 2B cluster energies. We expect MB cluster contributions 
containing or Vs^ij and Fg correlations to be important. In this work, as in earlier studies 
by FP and WFF, these contributions are calculated using chain summation methods [Q. 
The following three improvements were made by WFF in the basic VCS calculation used by 
FP: 

(i) The kinetic energy can be calculated using different expressions related by integra- 
tion by parts. If all MB contributions are calculated, these expressions yield the same 
result. However, they yield different results when only selected parts of the MB clusters are 
summed by VCS techniques. Studies of atomic helium liquids with VCS methods [^] find 
the exact result to be between the energies obtained using the Jackson-Feenberg (JF) and 
Pandharipande-Bethe (PB) expressions. Both the JF and PB energies were calculated by 
WFF, who took the average of the two as the result, and half the difference as an estimate 
of the error. 

(ii) The pair distribution function g'^{r) is defined such that pg'^ir) represents the prob- 
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ability of finding a particle at a distance r from a given particle. Conservation of number of 
particles then implies: 

h = p [dMl-g%r)) = l, (3.14) 



for both SNM and PNM. Also, since SNM has total isospin T=0, we have the following 
identity for expectation values of Tj ■ Tj operators: 

^- = t(0| E t.-t,|0) = -3. (3.15) 

The values of Ic and 1^ calculated using the VCS method are generally within a few percent 
of these exact results. However at small densities, the correlations become large due to the 
bound deuteron and virtual bound ^Sq states, causing Ic and I^ to deviate from their exact 
values by more than 10% in some regions of the dc, dt, a parameter space. Deviations of 
this size can also occur at large densities. Excursions into such regions of the parameter 
space are curtailed in SNM by minimizing: 



{H)+A 



{Ic - 1)^ + {\lr + 1)^' 



(3.16) 



with A chosen to keep Ic and Ir within 10% of their exact values. Only the A(/c — 1)^ 
constraint is applicable for PNM. 

(iii) In their calculation, WFF added the leading multiple-operator chain contributions 
to those summed via VCS. At po? these terms contribute less than 1 MeV/nucleon, while 
at 4po, they may contribute a few MeV/nucleon. Including these terms, WFF estimate the 
error in the calculation of the MB contributions due to kinetic energy, Vg and Fg to be ~ 
0.5 MeV/nucleon at po, and ~ 5 MeV/nucleon at 6po in SNM; the corresponding error in 
PNM is lower still. We note that these errors are negligible compared to the error implicit 
in the use of the simple given by Eq. (3.1). 

The computer program used to carry out the present calculations is built upon that of 
WFF and retains all the improvements made by them. 

The contribution of MB clusters involving spin-orbit correlations and interactions is 
calculated using methods developed by Lagaris Separable three-body clusters, with 

12 



correlations between pairs ij and ik, but not between jk, make the dominant contribution 
to (f + T)-MB-b, via v-ij or Vf terms in the Hamihonian. Chain diagrams with correlations 
between all three pairs ij,ik and jk were found to make a smaller contribution. We therefore 
sum all three-body separable contributions to {v + T)-MB-b, and estimate only the leading 
chain contributions. 

Lagaris also calculated part of the three-body separable contribution to (f )-MB-q using 
the U14 Hamiltonian The and (L ■ S)^ interactions that contribute to (f)-MB-q are 



significantly stronger in the A18 model than in the U14 model. We therefore include a more 
complete calculation of three-body separable terms and leading central chain contributions 
to (f )-MB-q in the present work. These calculations are outlined in appendix A. The present 
calculation of (f )-MB-q and the earlier calculation by Lagaris yield similar results for the 
U14 Hamiltonian: 1.61 (3.91) MeV versus 1.35 (3.0) MeV at kp = 1.33 (1.6) fm'^. The 
difference between the the two calculations is significant for the the A18+IX Hamiltonian, 
which makes much greater contributions to (f )-MB-q (Table |). 

We have summed the contributions of leading three-body and more than three-body 
diagrams to the expectation values of the V'^'" and static three-body interactions, using 
VCS methods described in Reference PUI. The A18 interaction induces stronger correlations 



than U14, as is evident by comparing T-2B-S for the two interactions in Table |. The larger 
correlations lead to comparatively larger {V'^'") and (^'^)- The {V^'^) in U14+VII and 
A18+IX models equal —3.49 and —3.60 MeV respectively at kp=1.33 fm~^, despite the 
smaller strength of V^"" in Model IX (^2^ = -0.0293 MeV) compared to Model VII (^2^ = 
-0.0333 MeV). The (v^) is larger for A18+IX (6.33 MeV at A;i.=1.33 fm'^) compared to 
U14+VII (3.99 MeV at A;f=1.33 fm^^). Only about half the increase in (v^) is due to the 
larger strength of in Model IX (f/o=0.0048 MeV) compared to Model VII (f/o=0.0038 
MeV). 

We note that the two-body correlations in this variational calculation do not have the op- 
timal form. The Fij are obtained by minimizing the sum of 2B contributions of the potential 
V (Eq. 3.3), with healing constraints imposed at dc and dt. More general correlations can be 
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generated by separately minimizing the two-body cluster contribution to each partial wave, 
specified by l,S,J and the relative momentum k Thus, the correlations f{l,S,J,k) 



depend on all the quantum numbers, and yield a lower 2B energy than the Fij operators 
defined in Eq. (3.2) for v with the same dc and dt. The MB contributions cannot be easily 
calculated with the general /(/, S, J, k), however. 

The Fij operators provide a good approximation to the two-body correlations in matter. 
The small differences between optimum /(/, S, J, k) and Fij can be accounted for by inclusion 
of the second order two-particle, two-hole contribution, AE2, calculated in correlated basis 
perturbation theory ||4l|j4^ . We estimate this contribution as described in Appendix B, 
approximating the AE2 by the difference SE2B between the 2B cluster energies calculated 
using /(/, S, J, k) and F^j. The values of a, dc, dt are determined by minimizing the energy 
calculated from the F^j, and SE2B, calculated for these optimum a,dc,dt, is perturbatively 
added to the energy. The calculated values of 6E2B (Table |) for the U14 Hamiltonian are 



slightly smaller in magnitude than the AE2 values reported in Reference for the same 
interaction. The 6E2B correction is larger for the A18 and A18-I-IX Hamiltonians, which 
predict stronger correlations in matter. The AE2 is known to be relatively larger for the 



A14 interaction E^], which induces stronger correlations than the U14. Our best estimate 



of the variational energy, given by ^{EpB + Ejp) + SE2B, is also listed in Table |. 

Over a decade ago. Day and Wiringa |^ calculated the ground state energies of SNM by 
means of the Brueckner-Bethe method. Contributions of up to four hole-line diagrams were 
included, in order to reduce the estimated error in the calculated energy to ~ ±0.18(1.3) MeV 
at kf = 1.33(1.6) fm~^. The calculations demonstrated that the three realistic, two-body 



potentials used in that work, Bonn (1975) |^3|,^, Paris |^ and A14, yielded essentially 
the same energy of ~ —14.9 MeV at kf = 1.33 fm^^. At kf = 1.6 fm^^, the Paris and 
Argonne V14 models yielded similar energies of ~ — 17.8MeV, while the Bonn result was 
somewhat higher, at ~ — 16.8MeV. Our results with U14 and A18 (Table |I|) are close to the 
earlier results at kf = 1.33 fm~^, and are about 1 MeV lower at kf = 1.6 fm"^. As with 
the results of Day and Wiringa, the present results show a remarkable model independence 
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of the energy of SNM for kf < 1.6 fm~^, calculated from realistic two-nucleon interaction 
models. The differences in the energies obtained using the different models are smaller than 
the estimated error in the many-body calculations. 

The energies of SNM and PNM calculated using the A18+IX Hamiltonian appear in 
Fig. |I] and Fig. along with the results obtained by WFF using A14+VII. As discussed in 
detail below, there appears to be a phase transition in both SNM and PNM with the present 
Hamiltonian. The curves marked LDP and HDP show the energies obtained for the low and 
high densitiy phases, respectively. 

At larger densities, the A18+IX energies are significantly higher than the A14+VII 
energies. A major part of this difference can be attributed to the difference between Vijk in 
Models VII and IX. For example, the contribution of Model IX exceeds that of Model VII 
by 8.2 (4.5) MeV to SNM (PNM) energy at p = 0.32 fm'^. The Model VII is unreahstic, 
however, as it overbinds '^H and ^He. The remainder of the difference between A18+IX 
and A14+VII is due to the stronger momentum-dependent part, Vg^ij of A18, which makes 
significant contributions at larger densities. 

In the early 1980s, FP added phenomenological density-dependent terms to the U14 
Hamiltonian, and adjusted their parameters to reproduce the empirical equilibrium density, 
energy and compressibility of SNM. These results also appear in Fig. |] and Fig. |^. Energy 
density functionals based on the FP results for E(p) of SNM and PNM reproduce the binding 
energies of nuclei from ^^O to ^°^Pb rather well [^]. 

The minimum energy calculated for SNM with the present variational wavefunction and 
the A18+IX Hamiltonian is —12 MeV, compared to the empirical value of —16 MeV. As 
previously mentioned, results of variational and exact calculations of light nuclei suggest that 
including three-body correlations in the wavefunction could lower the variational bound on 
the energy by more than 1 MeV, and that the true energy may be more than 1 MeV lower 
still. Thus, the underbinding of SNM due solely to deficiencies in the A18+IX model of 
nuclear forces is probably much less than 2 MeV. This underbinding is a very small fraction 
of the total potential energy of SNM at po = 0.16 fm~^, which is ~ —50 MeV for this 
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Hamiltonian. 

At approximately twice the equilibrium density, we observe what appears to be a tran- 
sition in SNM due to a change in dt, the range of tensor correlations. The energy of SNM is 
shown in Fig. ^ as a function of density and of dt/ro, where tq is the unit radius defined by 

yroV = l. (3.17) 

The energies in Fig. ^ have been minimized with respect to variations in the other two 
parameters, a and dc, at each p and dt/ro- For p < 0.32 fm~^, the minimum occurs at 
dt/ro ~ 4, whereas for p > 0.32 fm~'^, it shifts to dt/ro ~ 6. The upper line in Fig. ||, labeled 
LDP represents the minimum energies for dt/vQ ~ 4, while the lower line labeled HDP gives 
the energies of the minima at larger values of dt/ro- 

This transition is probably related to pion condensation [^,|1^,0, and its spin-isospin 



structure is discussed in the following sections. The transition does not occur in similar 
calculations of SNM using either the U14 or A14 and is thus sensitive to the two-nucleon 
interaction model. The Fujita-Miyazawa two-pion exchange three-nucleon interaction is a 
necessary ingredient for the transition to occur in SNM. A similar transition occurs for PNM 
(Fig. at a lower density (~ 0.2 fm~^) with both the A14+VII and A18+IX Hamiltonians, 
but it is not observed for U14+VII. The E{p, dt/ro) for PNM with A18+UIX appear in 
Fig. ^. Unlike the A14 and U14 interactions, the transition persists in PNM with A18, in 
the absence of the Fujita-Miyazawa interaction, although it occurs at a much higher density 
(~ 3.5po) in that case. 

The contributions to E{p, dt/ro) of SNM at selected values of p and dt/ro are listed 
in Table || for the A18+IX model. The interactions f*"^, v'^'^ and v'^, associated with the 
operators SijTi-Tj, cTi- a-jTi-Tj and 1, make the largest contributions to the energy of SNM. 
These contributions are listed separately, along with the contribution of the entire v. The v^'^ 
and v'^^ contributions come mainly from the one-pion exchange interaction. The expectation 
values of v'^'^ , V'^'^ and many-body cluster contributions are significantly different in the LDP 
and HDP. The pion-exchange contributions to the energies of PNM and SNM are listed in 
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Table |ITT| for the LDP and HDP. The OPEP used to calculate {v'^) includes a ttA^A^ dipole 
form factor with a 5 fm~^ cutoff as described in Section VI. The V'^'^ is that given by Model 
IX. The results indicate that for SNM the change in the pion-exchange contribution between 
the two phases comes overwhelmingly from the V"^"^ , which is required to produce the phase 
transition. However, for PNM a significant part of the enhancement in the pion-exchange 
contribution comes from v'^ , thus indicating the diminished importance of the three-body 
interaction to the phase transition, and the persistence of the transition in the absence of 

In the present calculation, the transition in SNM is of first order, as is evident in Fig. |l|. 
Thus, it is difficult to obtain information about the HDP by studying the LDP. The order 
of the corresponding transition in PNM is not as evident in Fig. |^, though it appears to be 
of first order as well. 



IV. PAIR DISTRIBUTION FUNCTIONS 

The two-body density f^{r), associated with the operator Of^-, is defined such that: 

(0| E B{n,)O^Q)=A( dhB{r)ff,{r), (4.1) 

for any function B(rjj). This relationship is used to calculate the expectation values of the 
pair interactions v^{rij)Ofj in the ground states of nuclei and nuclear matter, denoted by 
|0). The P2 associated with the static operators 1, Tj ■ Tj, CTi ■ crj, (Ti ■ (TjTi ■ Tj, Sij and 
SijTi ■ Tj, are denoted by pi, P2, p'2 1 p\ ^^^d p^^ ■ ^ properties of SNM discussed in this 
and the following two sections are obtained from these six two-body densities. In the case of 
PNM, we disregard isospin and consider only p\, p\ and p\. The central two-body density, 
P2, is proportional to the probability of finding a pair of particles separated by a distance r, 
and asymptotically approaches the matter density, p, as r^ 00. All other p2 ^^e associated 
with spin-isospin correlations and must therefore vanish as r^ 00. 

The static p\, calculated for SNM at the equilibrium density of p = 0.16 fm~^, appear in 
Fig. ^. The large magnitudes of the pg^^ indicate that short-range correlations in SNM are 
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strongly spin-isospin dependent. The nature of these correlations is more evident in the pair 
distribution functions Pr5A/('^)' which are proportional to the probability of finding a pair 
of nucleons with total isospin T, spin S and spin projection M, as a function of r. Forest et 
al. recently studied these densities in light nuclei |]^. 

The Pto 0(^)5 S=0,M=0 states, are spherically symmetric and are obtained from the 
P2's using spin-isospin projection operators, such as ^5=0 = |(1 — cTj • CTj)- In SNM, 

pfUr) = ^ m{r) + Pl{r) - 3p^,ir) - p^^(r)] , (4.2) 

pS,o(0 = ^ [P2(r) - Plir) - Pl{r) + pY{r)] . (4.3) 

Whereas in PNM, 

P'il{r) = \[pl{r)-pl{r)]. (4.4) 

The pair distribution functions in S=l states have a quadropolar deformation due to the 
tensor interaction. The pg m PNM are given by: 

pg(r) = Co(r) - 2C2{r)P2{cose), (4.5) 
p('ii(r) = Co(r) + C2{r)P2{cos9). (4.6) 

The expectation value of Bi{rij)Ps=i, where ^5=1 = ^(3 + cTi ■ a-j), can be calculated from 
either the two-body densities P2'-' or from the distributions pfj^- Equating the two results, 
we find: 

^o(r) = ^xl[3p^(r)+p^(r)]. (4.7) 

A similar calculation of the expectation value of BfSij yields: 

C2{r) = ^p^(r). (4.8) 

In the case of SNM, the Co{r) and C2{r) in T=0,1 states are obtained from the expec- 
tation values of Bi{rij)Ps=i and BtSij, multiplied by isospin-projection operators, Pr=o,i- 
For the T=0,S=1 distribution functions, we obtain: 



Coir) = ^ X 1 [Sp'iir) - 3p^(r) + p^(r) - pY{r)] , (4.9) 



C^ir) = — x- 
' 12 4 



P2W-PrW ■ (4.10) 



The corresponding T=l expressions are given by: 

Co{r) = ^ X 1 [9p^(r) + 3p^(r) + 3p^(r) + pr(r)] , (4.11) 



CJr) = -x — 
' 3 16 



3p2M+PrW • (4.12) 



The extrema of Ptim(^>^) fixed r occur for M = at 6' = and ^ = ^. These 



extrema are plotted as a function of r in Figs. and y. From Eqs. [4.5[ - [4.6| , it follows that: 

pPx±i{r, ^ = 0) = p5'io(^, ^ = ^/2), (4.13) 
P?,Ui(r, ^ = vr/2) = i (pg,o(r, ^ = 0) + p^^^ ,o(^, ^ = ^/2)) . (4.14) 

We also note that in SNM: 

(2) 2T + 1 

Pt!s,m(^ oo) = p . (4.15) 

The Po^] o(^) ^ = 0) ^/2) are shown in Fig. ^ for SNM at p = 0.16 fm^'^, along with the 
corresponding functions for ^H, "^He and ^^O, from Ref. 0]. At small r, this density is large 
for 6 = 7r/2, where the OPE tensor potential is attractive, and small for 6 = 0, where the 
OPEP is repulsive. Thus, according to Ref. [0, equidensity surfaces having Po^| o 

> 0.01 

fm~^ are toroidal in shape. The ratio Po'lo{r,6 = 7r/2)/pQ^| Q(r, 6* = 0) is a measure of the 
strength of tensor correlations in T=0 states. For the maximum possible tensor correlations, 
Po'}fi{r,6 = 0) is negligible compared to Po\,o{r-,0 = 7r/2). Fig. | thus indicates that the 
tensor correlations in T=0 states, in nuclei and in nuclear matter, have near the maximum 
possible strength at r ~ 1 fm. The peak value of Po^i^]w'('') is almost 2.5 times the asymptotic 
value of 0.01 fm"'^ at p = 0.16 fm~^. The spherically symmetric two-body densities for SNM 
and nuclei in the T=l, S=0 channel (pi^oo('')) shown in Fig. ^ These distributions peak 
at r ~ 1, where the nuclear force is most attractive, and the peak value is about 1.5 times 
the asymptotic value. Both Pi'^^ir) and pg^j j^.j{r) are supressed near r ~ by the repulsive 
core in the NN interaction. 
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According to Ref. 0, the Po^imI^") Pi^o,o(^) have universal shapes in light nuclei 
at small r. In Figs. ^ and we have scaled the densities in light nuclei such that their 
maximum values equal those of SNM. The two-body density distributions in SNM appear 
to have nearly the same shape as those in light nuclei for r ^ 1.5 fm. However, significant 
differences occur for r ^ 2 fm. Note that the average interparticle distance in SNM at 
p = 0.16 fm~'^ is also ~ 2 fm. 

Forest et al. argue in Ref. that the ratio, -R^d, of the maximum values of Po i m 
a nucleus A and the deuteron provides a good approximation to the Bethe-Levinger factor. 
La, of the nucleus A. The pion and photon absorbtion cross sections in light nuclei scale 
with Raci- The calculated value of Raci for SNM at equilibrium density is 1.59A, which 
corresponds to La =6.36. 

We observe interesting changes in the two-body densities between the low- and high- 
density phases in SNM and PNM. The existence of a pion condensate is indicated in the 
HDP, as discussed in the following sections. The Pr's m{^) SNM at p = 0.36 fm~^ in the 
LDP, with dt/vQ = 4, and in the HDP, with dt/ro = 6 are shown in Fig. p. The analogous 
p^]y^(r) in PNM at p = 0.20fm~^ are shown in Fig. ^. The differences between pair densities 
in LDP and HDP are more pronounced in PNM. In both PNM and SNM, we find that the 

(2) 

long-range part of the tensor correlations is enhanced. In PNM, the S=0 Pq q is suppressed 
in the HDP, whereas in SNM, Pi^o,o pg^] ^ are suppressed and Po^o.o is enhanced. 

In Migdal's approach |18|, the transition to the pion condensed phase is inhibited by a 
positive, short-range cTj ■ cr^Tj ■ Tj two-nucleon interaction f'^^(r), represented by the Landau 
parameter g'. The eigenvalues of cTj ■ cr^Tj ■ Tj are -3, -3, 1 and 9 in T,S=1,0; 0,1; 1,1; and 
0,0 states, respectively. Therefore, a strong positive f'^'^(r) favors the LDP, which has larger 
pair densities in T,S=1,0 and 0,1 states, and a smaller pair density in T,S=0,0. Similarly, 
in PNM, a large positive v'^, given by Eq. |2.13| , favors the LDP. The v'^'^{r) and v'^{r) in the 
U14, A14 and A18 models are shown in Fig. [1^. The positive v'^'^ and v'^ of U14 prevent 
a transition to the HDP in both SNM and PNM. In A14, v'^ becomes negative at small r. 
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and PNM is thus predicted by WFF to undergo a transition, while in A18, both v'^'^ and v'^ 
change sign and thus cannot prevent the transition in either SNM or PNM. The Urbana- 
Argonne potentials have similar forms, but are fit to different data sets. Whereas the U14 
and A14 were fit to n-p phase shifts available in the late 1970s and early 1980s, respectively, 
the A18 is fit directly to the 1994 Nijmegen p-p and n-p scattering database. Also, a much 
better fit was achieved by A18 (x^ per datum = 1.09). Thus, it is likely that A18 provides a 
more accurate representation of the NN interaction within the Urbana-Argonne framework. 

The f '^^(r) predicted by vr— and p— exchange potentials is positive for point-like nucleons. 
However, these meson-exchange potentials also contain a negative 5— function term. It is 
possible that this term is broadened by the finite size of nucleons, and that the total f°"^(r) 
changes sign at small r. The one-meson exchange representation of the NN interaction may 
not be reliable at small r, however. 

The v^\r)Sij Ti ■ tj interaction contains the main part of the OPEP. At r > 1 fm this 
interaction is essentially identical in the U14, A14 and A18 models. At r < 1 fm it is 
strongest in A14 and weakest in U14. 



V. ISOVECTOR SPIN-LONGITUDINAL RESPONSE 



Migdal IjTSl calculated the isovector spin-longitudinal (IVSL) response of nuclear matter 



using effective interactions. The IVSL response is defined as: 

RL{q,uj) = Y.\{I\0L{q)m^6{uj -uo-uj), (5.1) 
Ol{cO= E o-i ■ qr; ■ t e^'i-S (5.2) 

i=l,A 

where t is a unit vector in isospin space, and |0) and |/) represent the ground and excited 
states of the system, with energies ujq and uj, respectively. SNM has zero isospin, and Rl 
is therefore independent of the direction of t. The operator OL(q) represents the coupling 
of an external pion field to NM. In the case of PNM, we take t to be in the z-direction, 
such that T • t = —1. Ol{(i) then represents the coupling of a 7r° field to PNM. Migdal 
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assumed that this response would be dominated by a spin-isospin sound mode, and that the 
occurence of transition would be indicated by the vanishing of the corresponding excitation 
energy. 

Calculation of the response of NM from realistic interactions is an extremely difficult task. 
However, it is well known that the sums and energy- weighted sums of response functions 
are related to the two-body densities. For the IVSL response, the sum and energy-weighted 
sum are defined as: 



Jo 

Aq^Wiiq) = / RLiq,uj)ujdLU, 
Jo 



(5.3) 
(5.4) 



thus removing the dependence on the number of particles. In SNM they are given by: 



SL{q) = 1 + 1 I PVirMqr) ~ ptirMqr) 



p=2,6 



d^r. 



with the D^^{r) tabulated in ref |^^. In the case of PNM we obtain: 

1 



SL{q) = 1 



p2ir)3oiqr) - p2ir)j2iqr) 



d'^r. 



(5.5) 
(5.6) 



(5.7) 



and 



Dlir) 



Diir) 



3 
4 r 



v''{r){l- jo{qr)) -v (r)j2(gr) 



v\r){2 + jo{qr)) - (^"(r) - 2d\r))j2iqr) 



(5.8) 
(5.9) 



We note that terms with even p do not occur in PNM, and the are defined as in Eq. ( p.l2 - 
2.15|) . The energy- weighted sum of Eq. (|5.6| ) contains only the contributions from the static 
parts of Vij. These represent the dominant contributions to the energy- weighted sums of 



responses to electromagnetic fields ||48|j49|| , and presumably to Wl as well. 

The calculated values of Silq) in SNM at p = 0.16 fm~'^ and p = 0.36 fm~^, and in 



PNM at p = 0.16 fm and p =0.2 fm , appear in Fig. O. At the higher densities, results 
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obtained for both the LDP and HDP are shown. At equihbrium density, the Siiq) exhibits 
a small enhancement in the g ~ 2 fm~^ region. Indications of an enhancement of the IVSL 



response have been observed in (p, n) reactions |Q . In the LDP, this enhancement grows 
slowly with density. We predict a much larger enhancement at g ~ 1.3 fm~^ in the HDP. 

When the response is dominated by a single collective mode, i.e. when only one of the 
states |/) largely contributes to the sum in Eq. (5.1), the energy of the collective state is 
given by WL{q)/SL{q)- As an example, the energy of Feynman phonons in atomic liquid ^He 
can be obtained from the Wl/Sl ratio. The spin- longitudinal response of nucleon matter 
probably has a large spread in energy; nevertheless, we can define a mean energy of the 
response as: 

EM ^ . (5.10, 



The resulting values appear in Fig. [T^ for the cases considered in Fig [11. 

We note that in the LDP, -El (9) is larger than q^/2m in SNM as well as in PNM, 
indicating that the nuclear interactions push the response to higher energies on average. At 
g ~ the IVSL response is almost entirely due to spin-isospin correlations, and E^i^q ~ 0) is 
therefore large. In Migdal's picture 0], the energy of the collective spin-isospin sound wave 
with g ~ 1.3 fm~'^ decreases with increasing matter density, and pion condensation occurs 
when the energy vanishes. In the LDP, where dt ~ 4ro, the present calculations show that 
Eilq) does not decrease with density at any g. However, the HDP has a lower Eiiq) than 
the LDP in the vicinity of g ~ 1.3 fm~^. It is likely that a part of the response at g ~ 1.3 
fm~^ shifts to lower energies, or softens, as the system moves from the LDP to the HDP. 

VI. PIONIC INTERACTIONS AND EXCESS 

If the HDP is in fact a "pion-condensed" phase, it must have associated with it an 
enhanced pion field and enhanced pion exchange interactions between the nucleons. The 
interaction between two nucleons, due to the exchange of a pion of momentum q, is given 
by: 
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^S(q) = ■ q'^J ■ ■ e'^'-^^^-^^^A'iq), (6.1) 

m^[m^ + q ) 

where A(g) is the pion- nucleoli form factor. The expectation value of the interaction is 
trivially related to the sum Sl^q) of the IVSL response: 

T (T.<M)) = - 2, h"(SLiq) - l)A^(g), (6.2) 

T (E^S(q)) = — ./t^I 2, ^g'(^^(g) - 1)^'(^)' (6-3) 

^ \^ / rnliml + g^) 2 

in SNM and PNM, respectively. The expectation value of the complete v'^^ is obtained by 
integrating these expressions over dq: 



The function ^(g): 



dqiiq). (6.4) 



= :A^7%^?'|(1 - SL{q))A\q) (6.5) 



gives the pion-exchange nucleon-nucleon interaction (OPEP) contribution as a function of q, 
the magnitude of the momentum of the exchanged pions in SNM; ^(g) for PNM is diminished 
by a factor of 3. The calculated values of ^(g) for the cases previously discussed appear in 



Fig. 13 for SNM and PNM, where we have used: 



A(^^) = (6-6) 

with A=5 fm~^ as an illustration. In the LDP, the attraction from the OPEP comes from 
a broad region around g ~ 3 fm"^, whereas in the HDP, the attraction is relatively more 
concentrated at g ~ 1.5 fm~^. Also, much of the repulsion due to low g pions in the LDP 



is absent in the HDP. The total OPEP contribution, obtained from Eq. ( |6.4| ), is listed in 
Table |ITT[ In the case of SNM, the {v'") is not very different in the LDP and HDP; a much 
larger change occurs in the (V^^'^), also listed in table |ITT[ In PNM, however, the difference 
in (v^) is more pronounced between the LDP and the HDP. 
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The difference in the expectation values of the pion number operator in matter and for A 



isolated nucleons is called the pion excess |51| . The operator for excess pions of momentum 
q, due exclusively to OPE interactions, is given by: 



(6.7) 



The distribution of excess pions is given by: 



ri{q) 



sjml + q 



2' 



(6.8) 



which is shown in Fig. [T^for SNM and PNM. The 77(g) exhibits a sharp enhancement in the 
HDP at g ~ 1.5 fm^^. 

The integral of ^rS^^q)-, denoted by (^6n^^'^, is listed in Table |T|. The total {612^,) includes 
other contributions from the NN— i>NA and NN— A A transition interactions |^T|, which we 
have not calculated here. These additional contributions take into account the changes in 



the pion field due to NA- and AA-box diagrams and due to the V^^l 



VII. CONCLUSIONS 

We have studied the short-range spin-isospin structure of SNM and PNM, using one of 
the most accurate models of nuclear forces currently available. For SNM at equilibrium 
density, we predict a short-range structure very similar to that found by Forest et al. 
in light nuclei, at interparticle distances < 1.5 fm. Symmetric Nuclear Matter is bound 
by small localized regions of strong attraction in the NN potential in T,S = 0,1 and 1,0 
states. The two-nucleon densities are found to have large overshoots, ranging up to 2.5 (1.5) 
times the uncorrelated values (Figs. | and ^ in the T,S = 0,1 (1,0) attractive regions. The 
interaction, and consequently the pair density, in T,S = 0,1 states is highly anisotropic due 
to the presence of the OPE tensor force. The two-body density overshoots in this state 
have femtometer-sized toroidal structures similar to those found in light nuclei; we therefore 
expect this feature to occur in all nuclei. These short-range structures are not very sensitive 
to the uncertainties in models of nuclear forces as discussed in Reference 0. 
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We also find that the Argonne Vig plus Urbana IX Vijk model of nuclear forces, which 
offers one of the best fits to the Nijmegen NN scattering database, as well as to the binding 
energies of light nuclei |^ , predicts that both SNM and PNM will undergo transitions to 
phases with pion condensation at densities of ~ 0.32 fm~^ and ~ 0.2 fm~^, respectively. 
The occurrence of this transition is sensitive to the short-range parts of the CTj ■ ctj and 
cTj ■ (TjTi ■ Tj NN interaction, as predicted by Migdal The transition does not occur 
with the older Urbana and Argonne f 14 models in SNM, while in PNM it occurs with the 
Argonne v^, though not with the Urbana vu. 

It should be stressed that the present calculations of pion condensation in the HDP are 
incomplete. Although we have used a variational wavefunction of the same form to describe 
both the LDP and the HDP, a different form permitting long-range order should be used 
for the latter. A wavefunction describing a correlated liquid crystal, containing layers of 
spin-isospin ordered nucleons, may be more appropriate for the HDP. Takatsuka et al. 



have used such a wavefunction, denoted ASL for "alternating spin layers" . A correlated 
spin-ordered solid (SOS) wavefunction |B3| has also been used in the past. However, only 



lowest-order variational ||S3| or G-matrix [^] calculations have been possible with these 
wavefunctions having long-range order. In view of the small difference between the energies 
of the LDP, and the HDP, it may be advisable to use similar methods to calculate both. The 
chain summation methods used here cannot be used with either ASL or SOS wavefunctions. 
However, in future it may be possible to calculate energies of the LDP as well as correlated 
ASL and SOS phases, using cluster Monte Carlo methods |l6|,|35|,^. 



The indication of a phase transition obtained with the present calculation should be re- 
liable. The changes in the pair distributions, the IVSL response and the pion fields strongly 
suggest that the HDP will exhibit pion condensation upon inclusion of wavefunctions per- 
mitting long-range order. We hope that such a modification will not significantly alter our 
results for the short-range structure and energy of the HDP. The momentum of the con- 
densed pion field is ~ L4 fm^^ at the onset of the phase transition, which corresponds to a 
wavelenth of 4.5 fm. The interlayer spacing, which equals half the wavelength is thus 
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predicted to be ~ 2.25 fm in both SNM and PNM. This spacing is larger than the spacing 
of 1.5 (1.7) fm between the layers of a simple cubic solid at p = 0.32(0.2) fm~'^. 

Most calculations exploring the possibility of pion condensation in matter explicitly con- 
sider the ttNA coupling. The baryons in ASL and SOS matter are taken to be superpositions 
of nucleon and A-states. The NN — N A and NN AA transitions are considered in these 
approaches, which make the theory more difficult. In contrast, we do not consider the A- 
degree of freedom explicitly in the present work. The resulting effects, along with those due 
to other mesons and nucleon resonances, are implicit in the two- and three-nucleon poten- 
tials obtained by fitting experimental data. Our Hamiltonian has only nucleon degrees of 
freedom, and its predictions can be calculated using a variety of many-body techniques. As 
previously mentioned, however, the predictions regarding pion condensations are sensitive 
to the details of the short-range NN interaction. 

It is encouraging to note that our calculated density of ~ 0.32 fm~^ for the onset of pion 
condensation in SNM is within the range 0.32 to 0.48 fm~'^ favored by calculations using the 
ASL wavefunction |Q. As with WFF, we predict a lower transition density (~ 0.2 fm^^) for 
PNM. This is due to the fact that in most realistic models, v^{r), the relevant interaction 
in PNM, is softer than f°"^(r) (Fig. [TD|). Calculations with ASL wavefunctions predict a 
higher density of ~ 0.5 fm~^ for the onset of neutral pion condensation in PNM [^]. The 
possibility of charged tt~ and/or K~ condensation in neutron star matter at densities above 
our result of 0.2 fm~^ for 7r° condensation is currently being investigated by a number of 



researchers p5|-[57[l. 
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APPENDIX A: THE TREATMENT OF AND (L S)^ INTERACTIONS 

The leading contribution to (f )-MB-q is from the separable, direct three-body diagram 
shown in Fig. pTSja. The earlier calculation of (f )-MB-q by Lagaris |Tl| with the U14 potential 



included only the main term of this diagram, namely the term with only central correlations 
between interacting particles. While this approximation was justifiable for U14, the much 
stronger and (L ■ S)^ interactions of A18 (see Table |I[) require calculation of additional 
terms. In the present work we include all relevant terms of the direct diagram presented in 
Fig. pr5|a, and in addition consider: i) the interacting exchange three-body separable diagram 
(Fig. |TB|b), ii) the passive exchange three-body separable diagram (Fig. |TB|c), and iii)central 
chain diagrams with and without exchanges (Figs. |T5|d-g.) Terms in these diagrams are 
classified as either F-diagrams, in which the gradients in the interaction operate on passive 
correlations Fn, and K-diagrams, in which the gradients act on interacting correlations Fij 
and the uncorrelated wavefunction $. The direct three-body separable diagram ([l5|a) , has 
the general form: 

ki,kj,k; p,m,p' ,q,q' 

-(/^0^<0-40g)(/^0^/^'0^)]<I>3(^,j,/)rfV,,rfV,, (Al) 

with 

$3(Z,J,/) = ei(k,r,+k,.r,+k,.rO^ 

where C[ ] represents the so-called c-part, or the cr, r -independent part of the operator 
product, as described in Reference 0. Also, the Oij in the separated part of the expression do 
not operate on the fuOu. The expressions for the interacting and passive exchange diagrams 
(Figs. p!5| a and b) are obtained by replacing the uncorrelated wavefunction with 
either ^^{j, i, I) or i), and inserting the appropriate exchange operators (^ J2n=i,4: O^j 

or i '^ln=i,i^ii) to the far left of the operator product. 
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The three-body separable diagrams with m = 9 — 12 (Figs. |T5|a-c) have large F-parts, 



and smaller K-parts which we neglect. For the diagram presented in Fig. p!5| a we include 
terms with p, p' = 1 — 6 and QiQ' = 1 — 6. The dominant contribution comes from terms with 
p = 1 and g, q' representing either central or tensor correlations. The contribution from the 
diagrams presented in Figs. [T^ and |T3|c is somewhat smaller, and we only consider terms 
with p^p' = 1,2 and g, g' = 1 — 6. Contributions to the energy of SNM at p = 0.28 fm~^ for 



both the U14 and A18 models, without three-body interactions are listed in Table |V 



The (L ■ S)f,- in O™' ^^'^'^ can be decomposed as: 



(L ■ S)^_,. = -1(L ■ S),, + + ■ a,L' + ^a,,(L,L), (A3) 



where: 



aij{L, L) = ScTj ■ LcTj ■ L — cTj ■ (TjL^. (A4) 



The (L ■ S)ij part of v]^ ^^'^^ is treated along with the as described by Lagaris pT 



and includes all the diagrammatic terms therein. The parts are calculated along with 
the as described above. We also include selected aij diagrams, which are expected 

to make the most significant contribution to the energy. These contain direct F-diagrams 
(Fig. |l5|a) with: p'q'pq = cctc, tccc, ttct, cttt, tctc, and the exchange F-diagrams (Fig. p!5|b) 
with p'q'pq = cctc, tccc, where it is understood that 't' represents all tensor correlations with 
and without r,- - To or r,- ■ t/ factors. 



The central chain diagrams shown in Figs. p!5|d-g make a modest contribution to the 
energy. The K-diagrams are obtained by inserting appropriate chain functions into the two- 
body integrals. Three central chain functions, G^^,^^^ are defined in Reference 0, where 
xx' denote the exchange character of the two interacting vertices of the diagram. Each of 
these may be either direct (x = d), exchanged only with particles in the chain [x = e) or 
part of a circular exchange involving both interacting particles {xx' = cc). When considering 
only static interactions and correlations, the G^^, can be directly inserted into the two- 
body integrals. However, the gradients associated with momentum-dependent interactions 
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may operate on the Slater functions in the de and cc chains, thus yielding somewhat more 
complicated expressions. 

The direct central chain diagram (Fig. |T5|d) has neither interacting particle exchanged. 
The K-part of this diagram is obtained in the same manner as static diagrams, namely by 



dressing the two-body direct diagram with 



e'^dd - 1 



In order to calculate the de and cc diagrams (Figs, ^f and [T5|g), the operator product 
fijOfjV'^O'^fljOfj is written in powers of L^. The L° part does not contain gradients, and 
the associated de and ee chain diagrams are thus obtained as in the static case by dressing 
the direct two-body diagram with: 

(2GS, + G^^Gl + . (A5) 



The cc chain diagrams (Figs. |T5|e and |T5^) are calculated by replacing I'^j/s in the 
two-body exchange diagram by: 

le«- (/,, + sGl^f - S . (A6) 
s s 

Here, s is the degeneracy, and kj = likfTij) is the Slater function. The K-contribution of 
terms in the de and cc diagrams containing operators are calculated using the chain 
functions G^l and G^^ defined as: 

2 

g£ = / d\iF[^kiki (A7) 

2 

Gel = -^hj I dhiF^^kfUi (A8) 



with 

The chain functions result from the L'j- acting on the Fermi gas part of the wavef unctions. 
These depend on rj and due to the exchanges at the interacting vertices of the chain 
diagrams. 

The de contribution of the part is obtained by dressing the two-body direct diagram 
with: 
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(GS, + 2(5/A;JrJ)G£) , (AlO) 

while the cc contribution is calculated by replacing rijlijl[J s in the two-body exchange 
diagram by: 

,oi, [^iMl _ r^^^^Gl, + 2G^:) - . (All) 

The present calculation of (v)-MB-q includes terms involving only the static correlations, 
since they yield the dominant contribution to the energy. Without (L • S)ij correlations, 
the ffjOfjV^O^ffjOfj operator product does not contain terms with L"^^, and the above 
calculation of central chain K-diagrams is complete. 

APPENDIX B: PERTURBATIVE CORRECTIONS 

We wish to calculate the two-body cluster using a better variational wavefunction. The 
standard two-body variational wavefunction has the form: 

*2;. = ( E nr.j)0f,)<^2,, (Bl) 

p=l,8 

where $26 is a plane-wave Slater determinant. The are determined from the E-L equations 
obtained by minimizing the total two-body cluster energy 

C2 = ^ E(^2.|i/ - ^fcLl^2.) = ^ E E E - ^fcLl^2.), (B2) 

^m,n ^ km,kr, CTrr^CTn Tm,Tn '"' 

obtained after summing over k^, k„, a^, (Tn, Tm and t„ at chosen values of dc, dt and a. 
We propose instead to minimize C2{km,kn) separately in each partial wave, thus obtaining 
momentum- and channel-dependent correlations f{l,S, J, k). The perturbative correction to 
the energy in an /, 5*, J channel is then given by: 

SE^il, S,J)= E [C2{1, S, J; km, k^) - C2{1, S, J; km, (B3) 

where C2 is obtained from the f{l,S, J, k) and C2 is calculated with the operator f^O^. 
The uncorrelated two-body wavefunction, in relative coordinates, can be expanded as: 
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JOJ 



+ Ssi E \/47r(2J+l) i-^j,(fcr) :t^j{^(JlJM| JOIM) 

+ 551 E ^'-[ef_(r)3;jl,, + 0*l(r)3^i:i.], (B4) 



j>i 

7M 



where J± = J ± 1, J^^^j are spin-angle functions, {IS J M\lmiSms) are Clebsch- Gordon 
coefficients and: 



e2(r) = -V(2J+ + 1) (J+1JM|J+01M) i,^(/.r) . (B5) 
Tfie correlated two-body wavefunction can thus be written: 

*(r, A;, -5, M) = E V^SJ + l) i''[i?j,5=o(r) J^^oj + Rj,s=iir) yy^j{JlJM\J01M)] 

J 

+ E ^'-[<_(r)3;jl,, + Rl{r)yr^,j]. (B6) 



The wavefunction uj=i^s=o,i = '"-Rj.s in uncoupled channels satisfies the E-L equation: 



- —uj,s + [ + - \j,s{k)]uj,s = — A;'iij,5, (B7) 



and is normalized such that: 

uj,s{r^d,)^rj^{kr). (B8) 

The \j,s{k) are constants, which are varied to match the boundary conditions: 

u'j,s{^ = 4) = j,(/^r) + rj'^ikr), (B9) 

where o?^ = dc in the spin-singlet state, and dx — dt in the uncoupled spin-triplet state, 
which is affected by the tensor correlations. 

In the coupled channels, J = /±l,5' = l, we obtain a pair of coupled equations: 



u 



m m 2J-I-1 ' m 
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These equations have two sets of solutions, denoted by {u-,u!-) and {u+, u!+), with boundary 
conditions: 



u-{r = 4)= rjj_ (kr) ; a;_(r = dt) = 0, 

u'_{r = dx)= jj_ {kr) + rj'j_ {kr) ; uj'_{r = dt) = 0, 

where d^ = dc ior J = 1,1 = J_ = 0, ^Si state, and df otherwise, and 

M^(r = dt)=0 ; uj+{r = dt) = rjj^{kr), 

u'Jr ^dt)^0 ; u;'{r^ dt) = j, {kr) + rj' {kr). 



(B12) 



(B13) 



The A'^(A;), is adjusted to match the derivative boundary condition on the dominant wave 
{u- or CU+), while X^{k) is varied to match the zero derivative condition on the secondary 
wave {u+ and uj-). Thus, A'^(A;) and A* (A;) depend on J and the I of the dominant wave. 

The Rf_^ can be expressed as superpositions of the two solutions, which match the bound- 
ary condition: Rf^{r — > oo) = ^jj^{r — > oo). Evaluating the Clebsch- Gordon coefficients in 
ej^, we find: 

1 



M=±l 



R 



dM=±1 
Rj+ 



R 



M=0 



J^M=0 _ 



r 

1 

r 

1 

r 
1 



IJ + 1 



'J+1 



U-\/j + M+V J+1 

cu-y/j + UJ+VJ + 1 



(B14) 



With the wavefunction in hand, we can calculate the two-body cluster energy, and thereby 
determine 5E2{l,S,J). The differential equations for the wavef unctions are solved at several 
values of /c on a grid from to kf. The expressions for C2{1, S, J, k^, kn) are then integrated 
over the Fermi sea. In order to insure that the correlations functions, /, in coupled channels 
are positive, the solutions are matched to the asymptotic forms at the first node of the bessel 
function, rather than at the healing distance {dc or dt) for values of k large enough such that 
the node occurs within the healing distance. 
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The momentum-dependence of f{l,S,J,k) is not very large, however the dependence 
on /, S, J channels is significant. We find that no additional attraction is obtained in the 
^Si —^Di channel, and very little comes from the ^5*0 and other singlet channels as can be 
seen in Table |V[ This presumably indicates that the parameters dc-, dt and a are optimum 
for those channels. The bulk of the additional attraction due to channel- and momentum- 
dependent correlations is in the —^F2 and ^Pi channels, with somewhat less coming 
from the channel. The present f^O^ correlation operator, with d^ chosen according to 
Eqs. 3.8 and 3.9, is probably inadequate to simultaneously describe correlations in S- and 
P-waves accurately. However, the 6E2 correction has little effect on the critical density of 
the predicted phase transition. 
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FIGURES 

FIG. 1. E{p) of SNM calculated using A18+UIX. Included for comparison are previous cal- 
culations of E(p) using A14+UVII (WFF), and U-DDI (FP). The two sets of variational minima 
obtained at p > 0.28 fm"^ are labeled LDP and HDP. 

FIG. 2. E{p) of PNM calculated using AlB+UIX. Included for comparison are previous cal- 
culations of E{p) using A14-I-UVII (WFF), and U-DDI (FP). The two sets of variational minima 
obtained at p > 0.16 fm"^ are labeled LDP and HDP. 



FIG. 3. Constrained energies: Ec{p,dt/ro) = \{EpB + Ejp) + A [(4 - 1)^ + {\lr + if 



of 



SNM using A18-I-UIX. A was set to 1000 MeV, in order to keep the integrals of the two-body 
densities Ic and Ir within 5% of their exact values of 1 and —3 during minimization. 

FIG. 4. Constrained energies: Ec{p,dt/ro) = ^{Epb + Ejp) + A{Ic - if of PNM, using 
A18+UIX. A was set to 1000 MeV. 



FIG. 5. f^-^~^ for SNM at p = 0.16 fm-^ 

FIG. 6. Pt5m=o 1 o(''' ^ ~ 0)^/2) for SNM at p = 0.16 fm^'^, and for light nuclei scaled to 
match the maximum value of the SNM distribution. 

FIG. 7. Pt5m=ioo('") SNM at p = 0.16 fm~^, and for light nuclei scaled to match the 
maximum value of the SNM distribution. 

FIG. 8. Two-body densities Pt^s,m SNM at p = 0.36 fm^^. The full and dashed lines 
represent results for the LDP and HDP, respectively 

FIG. 9. Two-body densities pg'^j^ for PNM at p = 0.2 fm~^. The full and dashed lines represent 
results for the LDP and HDP, respectively 



FIG. 10. v°^'^{r) and v'^(r) in U14, A14 and A18 models of the nucleon-nucleon interaction. 



39 



FIG. 11. SL{q) for SNM and PNM. 



FIG. 12. WL{q)/SL{q) for SNM and PNM. q'^/2m (dash-dot curve) is plotted for comparison. 

FIG. 13. ^{q) for SNM and PNM. 

FIG. 14. Pion excess function, 77(g) for SNM and PNM. 

FIG. 15. (a)3-body separable direct, (b)3-body separable, interacting exchange, (c)3-body sep- 
arable, passive exchange, (d) direct central chain, (e) interacting exchange central chain (f) pair 
exchange central chain and (g) circular exchange central chain diagrams. 
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TABLES 

TABLE I. Contributions to the Energy of SNM in MeV 





U14 


U14 


A18 


A18 


A18-I-IX 


A18-I-IX 


A;/ (fm ^) 


1.33 


1.60 


1.33 


1.60 


1.33 


1.60 


p (fm ^) 


.16 


.28 


.16 


.28 


.16 


.28 


dc (rm) 


2.15 


1.79 


2.13 


2.08 


1.80 


1.67 


drt (fm) 


3.44 


2.86 


5.67 


4.76 


4.81 


3.81 


q; 


.80 


.80 


.65 


.61 


.80 


.90 


(T)-IB 


22.01 


31.85 


22.01 


31.85 


22.01 


31.85 


('u)-2B-s 


-56.40 


-75.21 


-66.17 


-100.68 


-66.41 


-105.34 


(T)-2B-s 


16.08 


21.72 


19.01 


28.51 


20.25 


33.47 


(w)-2B-b 


-3.23 


-6.37 


-2.29 


-5.02 


-2.38 


-5.76 


(r)-2B-b 


.80 


1.50 


.54 


1.10 


.62 


1.58 


(■u)-2B-q 


1.11 


3.56 


4.31 


11.39 


4.46 


12.44 


{v + r)-2B 


-41.64 


-54.80 


-44.61 


-64.71 


-43.46 


-63.61 


(■u + r)-MB-s 


4.76 


2.91 


6.47 


9.31 


5.50 


8.47 


\v + r)-MB-b 


-.21 


-1.02 


-.28 


-1.01 


-.22 


-.72 


(t;)-MB-q 


1.61 


3.91 


3.00 


8.54 


3.38 


10.62 




.00 


.00 


.00 


.00 


-3.60 


-9.96 




.00 


.00 


.00 


.00 


6.33 


22.09 


SE2B 


-0.95 


-2.06 


-1.30 


-2.35 


-1.89 


-4.07 


^{EpB + Ejf) + 6E2B 


-14.42 


-19.22 


-14.71 


-18.37 


-11.96 


-5.33 


^{EpB — Ejf) 


.30 


.27 


.55 


.82 


.60 


1.10 
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TABLE II. Contributions to E{p, df/ro) of SNM for A18+IX model in MeV 



p (fm 2) 


.16 


.28 


.28 


.36 


.36 


dt/ro 


4.20 


4.00 


6.00 


4.00 


6.00 


dc/dt 


.38 


.41 


.25 


.44 


.28 


a 


.80 


.90 


.96 


.98 


.88 


(T)-IB 


22.11 


32.10 


32.10 


37.96 


37.96 


{v + r)-2B 


-43.65 


-62.71 


-67.24 


-74.37 


-79.54 


{v + r)-MB 


8.71 


18.11 


35.29 


28.52 


43.70 


(T) 


42.26 


67.41 


72.28 


85.16 


85.97 


{v) 


-55.09 


-79.91 


-72.13 


-93.06 


-83.85 




-29.14 


-45.82 


-49.65 


-56.94 


-59.23 




-10.30 


-12.85 


-6.20 


-13.81 


-6.97 




-26.04 


-47.27 


-47.10 


-61.65 


-61.25 




-3.64 


-10.72 


-22.88 


-18.35 


-32.26 




6.42 


22.65 


23.76 


40.69 


41.71 


SE2B 


-1.91 


-4.65 


-2.26 


-6.78 


-6.20 


\{EpB + Ejf) + 


-11.96 


-5.22 


-1.22 


7.67 


5.37 


^(-E^PB — Ejf) 


.60 


1.05 


2.30 


1.45 


2.82 
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TABLE III. Expectation values of the Pion Exchange Interaction and Pion Excess Operators 



{v^)/A, MeV {V'^^)/A, MeV {6n);')/A 



SNM, p = 0.16 fm-3 


-31.53 


-3.64 


0.05 


SNM, p = 0.36 LDP 


-60.62 


-18.35 


0.09 


SNM, p = 0.36 fm-3, HDP 


-66.86 


-32.26 


0.15 


PNM, p = 0.16 


-0.67 


1.23 


-0.01 


PNM, p = 0.20 fm-3, LDP 


0.14 


1.90 


-0.01 


PNM, p = 0.20 fm-3, HDP 


-18.20 


-8.67 


0.05 


TABLE IV. MB-q Contributions 


in MeV to SNM at p 


= 0.28 fm~^ using 


U14 and A18 


NN-interactions 








diagram 


U14 




A18 


2B-dir 


4.18 




9.12 


2B-ex 


-0.62 




2.27 


(a) 


3.47 




8.57 


(b) 


0.22 
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TABLE V. Contributions to SE2 from l,S,J channels for A18+IX model in MeV 



p (fm 
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0.00 


-0.20 
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